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I review recent progress in understanding and modeling galaxy clustering in cosmological simu- 
lations, with emphasis on models based on high-resolution dissipationless simulations. During 
' the last decade, significant advances in our understanding of abundance and clustering of dark 

C^) \ matter halos allowed construction of accurate, quantitative models of galaxy clustering both 

\Q . in linear and non-linear regimes. Results of several recent studies show that dissipationless 

simulations with a simple, non-parametric model for the relation between halo circular vc- 
locity and luminosity of the galaxy they host, y max — L, predict the shape, amplitude, and 
luminosity dependence of the two-point correlation function in excellent agreement with the 
observed galaxy clustering in the SDSS data at z ~ and in the DEEP2 samples at z ~ 1 
over the entire probed range of projected separations. In particular, the small-scale upturn of 
' the correlation function from the power-law form in the SDSS and DEEP2 luminosity-selected 

samples is reproduced very well. At z ~ 3 — 5, predictions also match the observed shape 
q . and amplitude of the angular two-point correlation function of Lyman-break galaxies (LBGs) 

'_h ' on both large and small scales, including the theoretically predicted strong upturn at small 

\ scales. This suggests that, like galaxies in lower redshift samples, the LBGs are fair tracers 

of the overall halo population and that their luminosity is tightly correlated with the circular 
velocity (and hence mass) of their dark matter halos. 

1 Introduction 

During the last decade, large observational surveys of galaxies both at low and high redshifts 
have tremendously improved our knowledge of galaxy clustering, its evolution, and the relation 

between the galaxy and matter distributions. A coherent picture has emerged i n which bright 

galaxies are strongly biased with respect to the matter distribution at high redB hift^^^™™ 
and in which the bias decreases with time in such a way that the amplitude of galaxy clustering 
is only weakly evolving^ as expected in hierarchical structure formation 1 . 11 The bias is also 
in general scale-, luminosity-, and color-dependent. Brig ht (red ) galaxies are more stro ngly clus- 
tered than faint (blue) galaxies both in the local universd ^^^^ and in the distant pastC^EOlZI 

These trends are in general consistent with the picture in which galaxies reside in extended 
dark matter (DM) halos, forming via hierarchical collapse and merging of peaks in the ini- 
tial density field. Thus, for example, the stronger clustering of brighter galaxies can be read- 
ily understood if they tend to populate more massive halos, which are expected to be more 
clusteredPMSH Like galaxies, the halos are strongly clustered at high redshifts and their 
clustering strength evolves only weakly with tim d.^^ 

Clustering of halos of a given mass, formed in the standard structure formation scenario 
from a gaussian initial density field, is simple at large, linear scales, where it can be described 



by a single number — the linear bias@^ Recently, it has been shown that in addition to the 
mass dependence, the halo b ias depends on other halo properties, such as its formation time 
and mass concentrationP^^IEHI At small scales, the clustering of halos is more complicated 
and bias is scale-dependent, 11 the behavior resulting from dynamical evolution of halos in 
high-density environments?"!^ The non-linear clustering of halos on small scales and the pro- 
cesses shap ing it have been ext ensively investigated during the last decade using dissipationless 

simuiationsPiinEsinniiaiiia 

One of the most important theoretical advances of the last several years is development 
of the Halo Model (HM) framework, in which galaxy clustering on both linear and nonlinear 
scales is described quantitatively using spatial and mass distribution of DM halos calibrated 
against cosmological simulations and the Halo Occupation Distribution (HOD) — the probability 
distribution for a halo of mass M to host N galaxies with specified properties (e.g., luminosity, 
color, etc.). In its simplest form, the HO D is assumed to depend solely on the halo mass. This 
appears to be a fairly good approximation^^ However, in general the HOD can depend on other 
halo properties or large-scale halo environment. 

This framework proved to be extremely useful both in theoretical forecasts and interpretation 
of observed clustering data. For example, in the halo model the two-point correlation function is 
a sum of two separate contributions: the one-halo term, which arises from pairs of galaxies within 
the same dark matte r halo, and the two-halo term, which arises from pairs of galaxies from two 
different halos. 1 1 'The one-halo contribution dominates on small scales, while at scales larger 
than the size of the largest virialized regions clustering is due to the two-halo term. In general, 
the two terms are not expected to combine so as to give a perfect power-law correlation function. 
Departures from power-law are thus generically expected in this model. In addition, the model 
predicted that de viation s of the correlation function from a power-law should be even stronger 
at higher 

redshiftsPESl 

This is because at high redshifts the merger rate is higher and halos are 
more likely to have massive subhalos of comparable mass and luminosities. When the merger 
rate decreases, such massive companions disappear as they merge due to dynamical friction. 
These predicti ons have now b een convincingly confirmed by observations both at z = and 
higher redshiftsPMiranEHlEl 

Although the general idea of a galaxy-halo connection is definitely reasonable and is not 
disputed, the key question is how tight this connection is and whether properties of halos are 
tightly related to the properties of galaxies they host. Theoretical models of galaxy clustering, 
partially reviewed here, and their comparison to the wealth of current observational data start 
to shed light on these fundamental questions of galaxy formation theory. 



2 Theoretical Models of Galaxy Clustering 

Ultimately, one would like to simulate the distribution of galaxies in a large, representative vol- 
ume of the Universe, while reliably and self-consistently modeling their internal properties at 
the same time. Given that this is not yet feasible with the current state of our understanding of 
galaxy formation and the capabilities of the most powerful supercomputers, some phenomenolog- 
ical modeling and assumptions have to be made. Historically, galaxy clustering in cosmological 
simulations has been modeled using a variety of approaches. 

The most direct approach is to use cosmological simulations which include both dark matter 
and baryonic components, as well as galaxy formation processe s of radiativ e dissipation and 
phenomenological recipes for star formation and stellar feedback * ' ' ' ' ' Although proper- 
ties of galaxies are not yet modeled reliably in such simulations, they allow for unambiguous 
identification of galaxies as dense clumps of gas and stars, as well as measurement of basic galaxy 
properties such as stellar and baryonic mass, stellar ages, luminosities and colors. These observ- 
ables, in turn, allow for extensive comparisons with observations. The main disadvantange is 



Figure 1: Left: Dark matter distribution in a flat ACDM simulation of a 60ft -1 Mpc volume. Particles are 
color-coded on a greyscale according to the logarithm of the local density. One can see a network of filaments 
interconnecting groups and clusters. The network is filled with small, dense dark matter halos. Right: distribution 
of dark matter within the virial radius of a cluster-sized halo. The virial radius of the cluster is shown by the large 
circle. The volume enclosed by the cluster virial radius is filled with smaller dense subhalos, which are bound to 
the cluster and orbit within its potential. The circles indicate the individual objects identified by an automated 
halo finder (see Kravtsov et al. 2004 for details); the radii of the circles are proportional to the subhalo maximum 

circular velocity (r h oc V max oc M^ 3 ). 



that such simulations are generally computationally expensive, forcing one to sacrifice the size 
of the simulated volume or spatial resolution to make the simulations feasible. 

The most popular approach to modeling galaxy clustering employs semi-analytic modeling 1 
which uses phenomenological recipes for specifying when, where, and how galaxies form and 
evolve within dark matter halos, in conjunction with high- resolution dissipationless simulations 
modeling spatial distribution and merger histories of dark matter halos. Such hybrid models 
of increasing degree of sophistications have been used to model ev olution of ga laxy clustering, 
as well as trends with luminosity, color, and other galaxy propertiesP2ISlS1^212U These methods 
provide flexibility to explore the dependence of predictions on particular assumptions about 
galaxy formation physics, albeit at the expense of a fairly large number of free parameters, 
assumptions, and (often uncertain) parameterizations of the complex physical processes. 

The third, considerably simpler approach, is to use high-resolution dissipationless, DM- 
only simulations capable of following the evolution of both isolated, distinct halos and subha- 
los — the bound, self-gravitating dark matter clumps orbiting in the potential of their host 
hald . * I ' * * Subhalos are the descendants of halos accreted by a given system throughout 
its evolution, which retain their identity in the face of disruption processes such as tidal heating 
and dynamical friction (see Figure In the context of galaxy formation, there is little con- 
ceptual difference between halos and subhalos, because the latter have also been genuine halos 
and sites of galaxy formation in the past, before their accretion onto a larger halo. We thus 
expect that each subhalo of sufficiently large mass should host a luminous galaxy and this is 
indeed supported by self-consistent cosmological simulations^^ The observational counterparts 
of subhalos are then galaxies in clusters and groups or the satellites around individual galaxies. 
In this sense, I will use the term halos to refer to both distinct halos (i.e., halos not located 
within the virial radius of a larger system) and subhalos. 

With the assumptions that 1) there is one-to-one correspondence between halos in dissipa- 
tionless simulations and luminous galaxies and 2) there is a tight relation between halo properties 




Figure 2: Left: radial distribution of subhalos in high-resolution dissipationless simulations of cluster-sized 
halo. The solid line shows the density profile of dark matter. Different panels show radial number density 
profiles of subhalo samples selected using different criteria and different mass and circular velocity thresh- 
olds. In the top panels subhalos are selected using total bound mass and circular velocity at the present 
epoch, while in the bottom panels the mass and circular velocity values at the epoch at which each subhalo 
was accreted by the cluster are used. Right: radial distribution of subhalos in TV-body + hydro simulation 
of the same cluster with cooling and starformation. In this case, subhalos are selected using stellar mass of 
the galaxy they host. Comparison of the radial distribution in the left and right plots shows that at radii 
r/r v i r J> 0.15 — 0.20 the selection using circular velocity at accretion mimicks selection using galaxy stellar mass. 

Reproduced from Nagai & Kravtsov (2005). 



and properties of the galaxy they host, one can use halo distribution self-consistently modeled 
in such simulations to make detailed predictions for galaxy clustering. Confronting these pre- 
dictions with observations can test the validity of the assumptions above and constrain the 
relations between galaxies and their halos. In the remainder of this contribution, I will describe 
such studies in more detail. Their results show that this simple approach is remarkably accurate 
in describing the luminosity-dependence and evolution of galaxy clustering. 



3 Modeling galaxy clustering in dissipationless simulations 

3.1 Relating galaxies and halos 

Assuming that all luminous galaxies live in the dark matter halos identified in high-resolution 
dissipationless simulations, the key question is how we relate a galaxy of a given luminosity to 
a specific halo. In what follows, I will consider only luminosity among the possible properties 
of galaxies. Although it may be possible to relate other properties, such as color, to halo 
properties, these relations are considerably more uncertain and are likely to exhibit large scatter. 
Luminosities (especially in the red and infrared bands), on the other hand, are related to the 
total stellar mass of the galaxies. The stellar masses, in turn, can be expected to be related 
to the depth of the potential well of the halo and hence to its maximum circular velocity Knax 
(note that potential energy of a halo scales as W oc V^ ax ). This relation is thought to be the 
basis of the observed Tully-Fisher relation. Indeed, tight relations between halo circular velocity 
and stellar mass of the galaxies they host are expected in all models of galaxy formation. 

The halo property in simulations is attractive because it is unambiguous 

both for distinct halos and subhalos, which is not the case for the total mass. It should be 
noted that V mgx measured in dissipationless simulations will not correspond directly to observed 
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Figure 3: Bottom panel: The subhalo occupation distribution — the mean number of subhalos with maximum 
circular velocities V ma x > 70 km/s as a function of host halo mass. The solid line shows the mean total number 
of halos including the hosts (i.e., for each halo N is the number of subhalos plus one, the host halo itself), while 
the long-dashed line shows the mean number of satellite halos. Upper panel: the ratio of the square root of the 
second HOD moment to the first moment, a = (N(N - l)) 1/2 /(N), for the full HOD (solid points) and the HOD 
of satellite halos (open points). The dotted line at a = 1 shows the case of the Poisson distribution. Note that 
the HOD becomes sub-Poisson at small host masses. However, the HOD of satellites remains close to Poisson 
down to masses an order of magnitude smaller than for the full HOD. The dot-dashed line shows prediction for a 
of the total HOD, if the halo occupation distribution of satellite subhalos is described by the Poisson distribution. 

Reproduced from Kravtsov et al. (2004). 



rotation velocity of galaxies because dissipationless simulations do not take into account the 
effect of baryon condensation on circular velocity ?*"^" For our purposes, however, it is sufficient 
that a monotonic correlation between luminosity and halo Vmax is expected. 

The existence of such a monotonic relation allows for a simple, non-parametric model relating 
halo V max and galaxy luminosity. Specifically, the V max — L relation is derived from the measured 
abundance of halos as a function of their circular velocity and the requirement that the relation 
matches the observed luminosity function of galaxies. In the simplest version of the model, no 
scatter is assumed and the V max — L relation is derived by matching the circular velocity and 
luminosity functions: n(> Vmax) = n(> L). 

Such a simple model is bound to be too simplistic. First, there is certainly scatter between 
observed galaxy luminosity and their halo circular velocities. Such scatter can be added into 
the model at the expense of introducing a free parameter^! it is interesting that the amount of 
scatter between luminosity and Vmax may be constrained by joint comparisons of the model pre- 
dictions for galaxy-galaxy and galaxy-mass correlations because the two statistics have different 
sensitivity to scatter and this sensitivy depends on luminosity^! 

Second, circular velocity of halos measured at the epoch of observations is not expected to 
be the optimal choice for subhalos. For distinct isolated halos, the current circular velocity is 



a measure of their potential well assembled during evolution, and can therefore be expected to 
be tightly correlated with the stellar mass (or more generally the baryonic mass) of the galaxy 
the halo hosts. The circular velocity of subhalos in dissipationless simulations, on the other 
hand, is a product of both mass buildup during the period when the halo evolved in isolation 
and tidal mass loss, with an associated decrease of V m ^, after the halo s tarts to o rbit within 
the virialized region of a larger object and experience strong tidal forces* . * * * The stellar 
component of galaxies in centers of halos, which should be more tightly bound than halo dark 
matter, should be less affected by tidal forces and can stabilize the mass distribution (and hence 
Vmax) in the inner regions. We can therefore expect that luminosity and stellar mass of galaxies 
hosted by halos in dissipationless simulations should be correlated with the subhalo mass or 
circular velocity, V^ x , at the epoch of accretion, rather than with its current value. 

This i s b orne out by cosmological simulations that include gas dynamics, cooling, and star 
formation?!! Such simulations show that selection using V££^ results in subhalo distribution 
similar to the selection based on stellar mass of galaxies hosted by subhalos (see Figure [2J- One 
can therefore argue that a reasonable approach is to relate galaxy luminosity to the current halo 
circular velocity for distinct halos and to the circular velocity at accretion for su bhalos." The 
models based on the halo properties at accretion were recently used in several studie d. ^ I W I +V I W I O I 
The results discussed below show that such simple luminosity assignment model reproduces the 
luminosity-dependence of galaxy clustering at different epochs with remarkable, and perhaps 
surprising, accuracy. 



3.2 (Sub) Halo Occupation Distribution 

As I noted in § ^ the halo occupation distribution plays a central role in the modeling of 
galaxy clustering in the framework of the halo model. Given the assumption that every galaxy 
corresponds to a DM halo in high-resolution simulation, it is interesting to ask what form of the 
HOD the simulations predict. Analysis shows that for samples of halos with circular velocities 
larger than a c ertain threshold, the HOD is quite simple and can be parameterized with only a 
few parameter sJff^ l ^ 

Specifically, the HOD can be understood as a combination of the probability for a halo of 
mass M to host a central galaxy, associated with the halo itself and the probability to host 
a given number N s of satellite galaxies, associated with the subhalos (Figure 0). Such logical 
division makes physical sense, because central galaxies occupy a special location in the halo 
near the minimum of the potential well. Observational analogs of the central galaxies are, for 
example, the Milky Way with its system of satellites or a cD galaxy in a galaxy cluster. The 
HOD of the central galaxies can be approximated by a step-like function, while the satellite HOD 
can be well approximated by a Poisson distribution, fully specified by its first moment. The first 
moment of the satellite HOD can be well described by a simple power- law (N s ) oc M.P with (3 m 1 
for a wide range of number densities, redshifts, and different power spectrum normalizations. 

An important feature of the HOD shown in Figure |U is the "shoulder" near the minimum 
mass of the sample. The total HOD in the region of the shoulder is narrower than the Poisson 
distribution — - th e fact that produces a nearly power-law correlation function of galaxies at 
lower redshiftdA^ At higher redshifts, the shortening and steepening of the shoulder due to 
the younger age of host halos and more frequent presence of massive subhalos results in strong 
departures from power-law correlation function at small scales (see § \'AA\ and Figure EJ) • 

Remarkably, the form of the HOD derived for subhalos in dissipationless simulations is very 

"The distinction between circular velocity at accretion and current epoch has little effect at z > 1. This is 
because both the accretion and disruption rates are high at high redshifts. The accreted halos do not survive for 
a prolonged period of time, so that at each high-z epoch most of identified subhalos are recently accreted objects, 
which are yet to experience significant tidal mass loss. 




Figure 4: Left: comparison between the SDSS projected correlation function (points) and the correlation func- 
tion derived from halos (solid lines) for various luminosity threshold samples. For comparison we include 
the correlation function of dark matter particles (dotted lines) at the median redshift of the sample. Right: 
the meannumber of galaxies in halos of different masses for the luminosity samples shown in the left panels. 

Reproduced from Conroy et al. (2006). 



similar to galaxy HOD measured in iV-body+gasdynamics SPH simulations and models em- 
ploying semi-analytic models with dissipationless simulation jff l 65 l S 3 This supports the general 
framework of a close galaxy-halo connection. An additional important implication is that the 
physics shaping the HOD is relatively simple and is not sensitive to the details and specific 
assumptions of galaxy formation model. 

If halo circular velocity correlates with the luminosity of galaxies they harbor, we can expect 
that galaxy occupation distribution should have a similar form. Indeed, halo model fits to the 
galaxy clustering measurements in t he SP SS survey using the HOD form described above provide 
an excellent description of the dataP^E Moreover, the HOD paramete rs d erived from the data 
fits are in general agreement with the values expected from simulations^^ 



3.3 Modeling luminosity-dependent clustering in the SDSS survey 

Galaxy clustering as a function of luminosity at z = predicted using halo catalogs from 
dissipationless simulations and a simple, non-parametric model relating galaxy luminosity to halo 
Vmax) described in § 13.11 is shown in the right panels of Figure 0J The figure shows remarkably 
good agreement between predictions and clustering of the SDSS galaxies at all luminosities. Note 
that both the normalization and shape of the correlation function are reproduced. It is critical 
to realize that the agreement on scales r p ^ l/i -1 Mpc is due to the luminosity assignment 
scheme using V^ x . The luminosity assigned using present V max for subhalos would result in a 
significant under-prediction of amplitude of lo p at small scales, especially for fainter samples. 

At the same time, the figure shows that galaxies of different luminosity exhibit different 
amplitude and scale dependence of their bias with respect to the overall matter correlation 
function in simulation. Although the bias appears to be complicated and scale- and luminosity- 
dependent, it is faithfully reproduced by halos. 

The right panel of Figure 0] shows the first moment of the HOD (the mean number of 
galaxies in halos of mass M) for the luminosity samples shown in the left panels. The HOD 
shows explicitly that galaxies of higher luminosity reside in more massive halos. The overall 
shape of the HOD, however, is similar for all luminosities. 
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Figure 5: Left: comparison between the DEEP2 projected correlation function (points) and the correlation 
function derived from halos (solid lines) for various luminosity threshold samples. For comparison we include the 
correlation function of dark matter particles (dotted lines) at the median redshift of the sample. Right: similar 
comparison for the angular correlation function of the Lyman Break Galaxies with different apparent magnitude 
limits in the Subaru survey at z ~ 4. Note that the galaxy correlation function strongly deviates from the power 
law at small scales, the behavior expected for high- z objects. The dotted lines show the corresponding correlation 
function of matter in the simulation. At these high redshifts the LBG galaxies are highly biased with respect to 
the matter distribution. Reproduced from Conroy et al. (2006). 



3.4 Modeling galaxy clustering through cosmic time 

Figure |5] shows similar comparisons with galaxy clustering as a function of luminosity at z ~1 
in the DEEP2 survej 5 ^ and clustering of Lyman Break Galaxies (LBGs) in the Subaru survej^ 
at z ~ 4. At both redshifts, the overall agreement is again excellent on all scales. Small 
discrepancies at r p <^ 0.5/i -1 Mpc for the Mb — 51og/i < —20.5 DEEP2 sample may be attributed 
to cosmic variance and poisson noise. 

It is worth stressing again that this remarkable agreement between observed and model 
clustering is achieved using the halo distribution in dissipationless simulations with a simple, 
non-parametric relation between galaxy luminosity and halo circular velocity. The luminosity- 
dependent bias at all redshifts hence seems to be driven entirely by the fact that brighter galaxies 
reside in more massive halos, with the correspondence between halo and luminosity determined 
by matching the observed luminosity function to the dark matter halo circular velocity function. 
This may appear as a reasonable and not unexpected result for lower redshift galaxies where 
luminosity may be expected to be a good tracer of stellar mass. The agreement for the LBG 
galaxies is more surprising and suggests that, like galaxies in lower redshift samples, the LBGs 
are fair tracers of the overall halo population and that their restframe UV luminosity is tightly 
correlated with the circular velocity (and hence mass) of their dark matter halos. 

Note that at z ~ 4 the correlation function of the LBG galaxies exhibits strong departures 
from the power law at sm all scal es, the behavior that was expected based on the halo model 
arguments and simulations?^"^ Similar behavior was also recently observed in clustering of 
LBG galaxies in the GOODS survey at z = 4 and z = ty^ 1 

The analysis of the HOD of the z = 4 halos in the simulation supports the model in which 
most LBGs are the central galaxies in their host halos with luminosity tightly related to the 
halo circular velocity and mass. Most LBGs have no neighbors within the same halo. However, 
a fraction of them do and it is this fraction that is responsible for the strong upturn in the 
correlation function at small scales. The reason that small-scale upturn in the correlation func- 
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Figure 6: Comparison of the projected two-point correlation function for SDSS galaxies with luminosities M r < 
—20 and predictions from dissipationless simulations with different assumptions about fraction of orphan satellite 
galaxies (galaxies with subhalos disrupted by tides, See § 13.51 . The solid line black line consistent with the data 
is the prediction without orphans (i.e., using only subhalos in the simulation; this is the same line as in Figure0J. 
Other lines correspond to orphan fraction increase by 10% (the lines thus show effect of orphan fractions between 
10% and 50%). Orphan fractions of > 10% modify the amplitude and shape of w p sufficiently to break the 
agreement with the data. Figure courtesy of S.C. Conroy. 

tion becomes more pronounced at higher redshifts is that a larger fraction of the high-z halos 
has a satellite of comparable mass relative to the lower- z objects. In terms of the HOD, these 
differences manifest in a shorter and less flat "shoulder" of (N(M)) near the minimum mass of 
the sample at higher redshifts. 

By accurately reproducing both the small-scale upturn in oo{9) and the large-scale clustering, 
our model accurately predicts not only the correct distinct halos to associate with LBGs (the 
'2 halo term' in halo model jargon) but also the number of LBGs within a distinct halo (the 
corresponding '1 halo term'). 

3.5 Are we missing galaxies in dissipationless simulations? 

In dissipationless simulations, subhalos experience tidal heating and mass loss by the tidal field 
of the host halo. Such tidal mass loss and disru ption ar e not numerical artefacts but are real 
physical processes affecting evolution of subhalos 1 . 1 11 However, one can argue that baryonic 
dissipation effects may greatly enhance resistance of galaxies and their halos to tidal disruption. 
A decision thus has to be made as to subsequent evolution of galaxy when its subhalo is disrupted 
by host halo tides . b One extreme is to assume that stellar components of galaxies are always 
resistant to tidal disruption and are thus never disrupted. Such an assumption leads to up to 



6 Operationally, disruption should be understood as a significant tidal mass loss which brings the total mass of 
a halo below a threshold of a halo sample or resolution of the simulation. 



~ 30 — 40% of satellite galaxies without subhalos in dissipationless simulations . I I Another 
extreme is to assume that stellar systems get tidally disrupted at the same time or shortly after 
the host halo is disrupted (as is assumed in the model used for comparisons with data above). 

In order to assess the quantitative impact of such possible "orphan" galaxies on clustering 
statistics, the following simple test can be performed. The fraction of subhalos in a simulation is 
increased in such a way that the subhalo occupation function, (N(M)) sa _ t , increases in amplitude 
while maintaining the same shape. In order to match the observed galaxy luminosity function, 
the overall V max threshold for a given sample is simultaneously increased, such that the number 
density does not change. Figure El shows the predicted projected 2-point correlation function of 
galaxies with M r < —20 (results are similar for other luminosities) with different assumed frac- 
tions of orphan satellite galaxies. Even a small, ^ 10%, contribution of orphan galaxies breaks 
the excellent agreement of the model with the data. Note that the amplitude of the correlation 
function is more sensitive to the presence of orphan galaxies at small scales (r p <^ 1ft, -1 Mpc). 
Thus, addition of a significant fraction of the orphans changes not only the amplitude, but also 
the shape of the correlation function. 

The existence of a significant (~ 30 — 40%) fraction of orphan galaxies in the real universe 
would imply that the excellent agreement with the observed amplitude and shape of the corre- 
lation function at different luminosities and at different redshifts is fortuitous. The disruption 
of stellar systems by tides is therefore an important problem which should be tackled with 
high-resolution numerical simulations that include stellar component. 

3.6 Future prospects 

The presented results indicate that the clustering can be modeled quite successfully with dis- 
sipationless simulations. As the amount and quality of data on galaxy clustering continues to 
improve, we can expect more stringent model tests. It would not be surprising if data soon 
reveals limitations of such a simple model. As I noted above, an interesting avenue is to per- 
form joint comparison of the model to the galaxy-mass and galaxy-galaxy correlations. Such 
comparison can potentially constrain the level of scatter in the L — V max relation. 

On the theoretical side, comparisons with data can be done with simulations of larger vol- 
umes to reduce statistical and cosmic variance errorbars. Smaller error bars should allow us 
to see subtle deviations from the data. At the same time, the L — V max model used in dissi- 
pationless simulations can be directly and thoroughly tested with iV-body+hydro simulations 
with cooling and star formation. One important question that can be addressed in simulations 
is tidal disruption of galaxies and their DM halos and existence of "orphan" galaxies. This work 
is currently underway. 
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